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Abstract. We investigate the asymptotic behaviour of large eigenvalues for a 
class of finite difference self-adjoint operators with compact resolvent in I 2 . 



1. Introduction 

1.1. General remarks. Infinite tridiagonal tridiagonal matrices called "Jacobi matri- 
ces" have been investigated in many recent papers in relation with various questions of 
pure and applied mathematics (see [T1I MTuhT2] ) . In [Tll^IISllMTB] the authors investigate 
Jacobi matrices acting in I 2 as unbounded self-adjoint operators with discrete spectrum 
and asymptotic formulas for large eigenvalues are given. This type of analysis is of par- 
ticular interest in Quantum Physics when information about physical parameters can 
be deduced from the spectral asymptotics of a concrete model. 

Except for recent work [7] there are no corresponding work concerning the asymp- 
totic analysis of large eigenvalues of higher order symmetric difference operators. This 
fact is not so surprising because higher order difference operators had not been studied 
from the spectral point of view, up to the last few years. On the other hand there are 
already works dealing with spectral properties of difference operators of higher order, 
see, for example [3JIHH]. It is natural to ask if known results on tridiagonal matrices 
can be generalized to higher order difference operators and to look for applications. 
As an example of possible application let us mention here the problem of the behavior 
of large singular values for a non-symmetric Jacobi matrix J with discrete spectrum. 
Indeed, these singular values are eigenvalues of J* J which is a symmetric difference 
operator of order four. 

The aim of this paper is to obtain a simple remainder estimate in the asymptotics of 
eigenvalues for a large class of symmetric higher order difference operators. Applying 
this result to tridiagonal matrices, 

(i) we obtain asymptotic estimates of eigenvalues for Jacobi matrices which cannot 
be treated in [HHEHS] (see Section O]). 

(ii) but with remainder estimates less precise than those of [TJ[2j[5l|6] . 

Since our assumptions are weaker than in [TJEJUHI]) we must overcome some additional 
difficulties, but the main idea of our approach remains the same as in [2]. Although the 
approach of this paper is used to obtain the simplest remainder estimate, it is possible 
to follow the idea of [5] in order to compute further terms of the asymptotics with 
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smaller remainders under stronger "conditions of smoothness" imposed on the entries. 
It is an open question how to extend the methods of the present work to the case of 
non-smooth entries. Finally notice that we (and most of the above authors) have not 
studied asymptotic properties at infinity of the eigenvectors of A. 

1.2. Main statement. Let I 2 — l 2 (N*) be the Hilbert space of square summable 
complex valued sequences x : N* — > C equipped with the norm 

w[ P = (Ei a; 0')f' ~ 

We fix d: W ->• K and introduce 

OO 

V:=[xef: ^|d(j>(j)| 2 <oo}. (1.1) 
j=l 

Then we consider the self-adjoint operator D : T> — > I 2 given by the formula 

(Dx)(j) = d(j)x(j) for x e V (1.2) 
and a finite difference operator A' : V — > I 2 of the form 

(A'x)(j)= J2 (ai(j)x(j + l) + ai(j-l)x(j-l)) (1.3) 

l<Z<m 

where the coefficients a; : N* — >• M, I = 1, . . . , m satisfy 

I°'Cj--0] + I°«(j-)I , . (i.4) 

We assume ai(j — I) = = x(j — I) when j < I in (jl.3[) . We investigate the operator 

A = £> + 74' (1.5) 

under the additional assumption 

d{n) > oo, (1.6) 

n— >oo 

which ensures that ^4 has compact resolvent, hence there exists an orthonormal basis 
(w„) ne N* such that Av„ = X n (A)v n holds for n € N*, A n (A) — > oo as n — >• oo and 
(A n (A)) nS N* is arranged increasingly, i.e., X n (A) < X n+ i(A) for any n e N*. 

Assumption (HI). The off-diagonal entries ai(n), 1 < I < m satisfy the asymptotics 

ai(n) = cin s ' + O^" 5-1 ) as n -> oo, (1.7) 

where 5/, q, ? = 1, . . . , m are some fixed real numbers and S > max{<5i, . . . , <5 m }. 

Assumption (H2). The diagonal entries d(n) satisfy the asymptotics 

d(n) = c n 5a + en 5 "- 1 + 0(n S °- 2 ) as n -> oo, (1.8) 

where 5o > 0, Co > and c S K are fixed. 

Our main result is the following 

Theorem 1.1. Let A = D + A' be defined by (|l.i p -(|1.6 |) . If both assumptions (HI), 
(H2) hold and if k := 5q — S > 0, then 

X n (A) = d(n) + 0(n s - K ) as n -> oo. (1.9) 
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1.3. Comments. 

a. If (jTTTJ) is replaced by the weaker condition 

ai{n) = 0(n s ) 

then the min-max principle allows us (see Theorem 13. ip to prove 

X n (A) = d(n) + 0(n s ) as n -> oo. (1.10) 

The main purpose of Theorem 1 1.1 1 is to show that it is possible to replace the estimate 
(|1.10j) by the improved estimate (|1.9j) . 

b. For any fixed jgZ the assumptions of Theorem 11.11 imply 

ai(n + j) =Q(n -«) as „^^ (LU) 
d(n) 

We observe that the assertion of Theorem 1 1 . 1 1 holds for any fixed re > while all papers 
PMUHT] assume re > 1. 

c. We observe that 

d(n + 1) - d(n) ~ <5 c n' 5o_:L as rt ->■ oo (1-12) 

and we can treat the case < 6q < 1 when d(n + 1) — d(n) — > as n ~ > oo, while the 
papers [TJ[SHZ] assume that 

liminf (d(n + 1) - d(n)) > 0. 

d. Thcorcm ll.il will be obtained as a special case of more general estimates described 
in Section [5] (see Theorems 15.11 and 15.21) . In Section 15.41 we give asymptotic estimates 
of eigenvalues for some cases of not power-like entries. 

1.4. Contents. In Section [5] we check that the operator A is well defined under as- 
sumption (|1.4[) and its resolvent is compact under assumption (|1.6[) . 

In Section [3] we show how the min-max principle ensures the estimate (jl.lOj) if (jTTTJ) 
is replaced by the weaker condition ai(n) = 0(n s ). 

In Section [4] we present basic ingredients of our approach based on the construction 
of operators which are unitarily similar to A with smaller off-diagonal entries. A similar 
idea is often used to investigate spectral asymptotics of self-adjoint problems defined 
by a linear PDE, e.g., in relation with the semi-classical approximation in Quantum 
Mechanics. 

In Section ISTTI we state a "general estimate" (Theorem l5.1[) . In Section [5721 we derive 
Theorem 15.21 which is an application of the general estimate to power-like entries. In 
Section 15.31 we easily check that Theorem 11.11 is a special case of Theorem 15.21 and in 
Section 15.41 we apply we apply Theorem 15.11 to cases where the entries have different 
asymptotic behaviors. 

Finally, we complete the proof of the general estimate in Section [6] 

The main result (Theorem 1 1 . 1 [) stated above derives from Theorem 15.11 as follows: 

Theorem lS.ll l 



Lemma 14.2 



Theorem [572] =^> Theorem O 
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2. The operator A 

Proposition 2.1. Let A = A' + D be the operator defined by - ()1.3[) and (|1.5[) . 

(i) // (jl.4p holds, then A is self-adjoint. 

(ii) // moreover (jl.6D holds, then A has compact resolvent. 

It is well known (see [5]) that (|1.5JI defines a self-adjoint operator T> — > I 2 provided 
A' has zero relative bound with respect to D, i.e., if for any e > there is C e > such 
that 

\\A'x\\ P < e\\Dx\\ P + C e \\x\\ P for xeV. (2.1) 
Before starting the proof of Proposition 12 . 1 1 we introduce some notations. We recall 



that the scalar product in I 2 is defined by (x,y) — Y^kLi ^(k)y(k) an d we denote by 
( e n)^Li the canonical basis of I 2 , i.e., e n (j) — 6j >n where 5 n> „ = 1 and 5j. n = for 
3 n. 

Then we observe (see [8]) that it suffices to show (|2.1[) for x G c o, where c o is the 
linear subspace of I 2 generated by the canonical basis, i.e., 

x G c 00 <=> #{j : x(j) ^ 0} < oo. (2.2) 

We denote by B(Z 2 ) the algebra of bounded linear operators on I 2 with the norm 

||T||= sup \\Tx\\ l2 . (2.3) 
lkll)2<i 

The shift operator S G S(i 2 ) is defined by 

S'e„=e„ + i (2.4) 
and for any o:N*->Cwe denote by a(A) the closed operator in I 2 given by 

(a(A)x) (j) = a(j>(j) for x G c 00 . (2.5) 
We can then rewrite the definition of A' in the form 

A'x= (S l ai{A)+ ai (A)S l *)x (2.6) 

l<l<m 

where x G coo and is the adjoint of S l . 

Proof, (i) It suffices to show (12.1[) for x G coo- For ( = l,...,mwe denote 

A\ :=5 i a i (A)+a i (A)5 i *. (2.7) 

For arbitrary ir G Coo we can write 

\\A\xfp < (WS'a^xWp + \\ai(A)S* l x\\ P ) 2 . (2.8) 

Since S l ai(A)S* 1 — a;(A — /) with the convention that ai(j — I) = if j < I and S l is 
an isometry, the right-hand side of (|2.8|l can be estimated from above by 

2||5 i a i (A)a;|| 2 2 +2||a ; (A)^x|| 2 2 = 2| | ai (A)a;| | 2 2 + 2||a z (A - /)x|| 2 2 . (2.9) 

Then taking y = (i + D)x we obtain 

II^Ki + ^r^ll^EiMiMi)! 2 (2-10) 
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with 



For N e N* we denote by Hjy the orthogonal projection onto {e ra }i<„<jv and 11^ 



J-ILv. Then (|2~TTj) implies ||A'(i + I>) _:1 n^|| -» as TV -> oo. Thus for a given e > 



we can find iV(e) G N such that 

\\A' t (i + Dy'U'^yWp < sWvWp 

and we deduce ||AjJc||j2 < e||(£) + i)LT^/ e sa;||;2 + ||Ajn^( e )|| ||x||{3. 

(ii) The operator Aft + D)^ 1 is compact as limit of finite rank operators A[(i + 
D)- 1 ^ in the norm of B{1 2 ). Then (i + A)' 1 Aft + D)- 1 = (i + D)- 1 - (i + A)- 1 

is compact and compactness of (i + D) -1 (due to (|1.6p ) implies that (i + A)^ 1 is 

compact. □ 

3. ASYMPTOTICS BY MIN-MAX PRINCIPLE 

3.1. Statement. In what follows for a sequence x(n) we will use the notation 

(Ax)(n) := x(n + 1) - x{n). (3.1) 
The purpose of this section is to prove the following 

Theorem 3.1. Let A = D + A' be given by (ll.lj) - (|1.6l) . Assume moreover that there 
exist C > 0, (5 6 8, k > satisfying S + k > and no £ N smc/i £/ia£ 

(Ad)(n) > C" V+"- 1 /or n > n , (3.2) 

and, /or £ = 1, . . . , m, 

ai(n) — 0(n s ) as n — > oo. (3-3) 
TTien one /ias £/ie Zarge n asymptotic formula 

\ n (A) = d(n) + 0(n 5 ). (3.4) 
We observe that due to f|3.2[) there exists c > and ni such that 

d(n) > cn 5+K for n > n% (3.5) 
and with ((331) imply (iLTTT) . 

3.2. Auxiliary estimates. 

Lemma 3.2. Assume that a: N* — > R satisfying the two conditions 

a(j)> Yl (|oi(j)l + |ai(j-OI), (3-6) 

l<i<m 

^4 ~~ ^ as j — >■ oo. (3.7) 

T/ien /or every b£N* i/ie estimate 

d- < X n (A) < d+ (3.8) 
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holds with 

d-:=mt{d(j)-a(j)}, (3.9) 

]>n 

d+ :=sup{d(j)+a(j)}. (3.10) 

j<n 

If moreover there exists jo £ N* such that 

\(Aa)(j)\ < (Ad)(j) for j > j , (3.11) 
then there exists ri\ £ N* such that 

\X n (A) - d(n)\ < a(n) for n > n x . (3.12) 

Proof. Let V n denote the linear subspace generated by {£j}i<j< n anc ^ Ki" denote its 
orthogonal complement in I 2 . Then suitable versions of the min-max principle give 

inf (Ax,x) < X n (A) < sup (x,Ax). (3.13) 

\\x\\ 2 <l N| 2 <1 

Let A\ be as in (|2.7p . Then writing 

(x,A[x) = ^a l (j)x(j + l)xJj)+^a l {k- l)x(k - l)~xjkj (3.14) 

3 k 

with k = j + I we can estimate | (a;, A[x) | by 

X)2|a,0-)l Hi + l)x(j)\ < E MMWJ + 0I 2 + k(j)l 2 )- (3-15) 

Therefore the right-hand side of (|3.15p can be written in the form 

E + oi 2 + E |aiCj)ii*Cj)i 2 = E(i°'0' - oi + m)\mj)\ 2 (^) 

k j j 

and we obtain 

\(x,A'x)\ <E«0')Ki)| 2 - (3-17) 

i 

Next we observe that (|3 . 1 T[) implies 

(x, (d(A) - a(A))x) < (x, Ax) < (x, (d(A) + a(A))x) (3.18) 
and using (|3.18p we can estimate the right-hand side of (I3.13[) by 

sup (x, (d(A) + a(A))> = d+ 

xEV n 

M 2 <i 

To complete the proof of (|3.8[) note that 

inf (x, (d(A) - a(A))x) = d~ 
NI 2 <1 

is smaller than the left-hand side of (|3.13p . In order to show (|3.12p we observe that 
(|3.1ip implies 

d(j) - a(j) < d(j + 1) - a(j + 1) for j > n , (3.19) 
d(j) + a(j) < d(j + 1) + a(j + 1) for j > n , (3.20) 
and consequently d^ = d(n) ± a(n) for n. > m, hence (|3.12[) follows from (|3.8p . □ 
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3.3. Proof of Theorem [1ITJ 

Proof. Due to (TETu]) there exists C > such that {jDjJl , (|X71) hold with 

a(j) := C / (3.21) 

and (11.9[) ensures the estimate 

|(Aa)(j)| ~ I^Co/- 1 < \S\ C Cr K (Ad)(j) for j > j . (3.22) 

Since |<5|CoCj~ K — > as j — > oo, it is clear that (|3.22p implies p. lip if jo is large 
enough. Thus p. 121) holds with a(j) given by (|3.21[) and the proof of (|1.12p is complete. 

□ 

4. Basic ingredients of the approach 
4.1. Main ideas. We write the following formal development of the conjugate 

B n := e- iP "Ae ip " = A + [A,iP n ] + ±[[A,iP n ],iP n ] + ... (4.1) 

where P n is self-adjoint and of finite rank for simplicity. Then X n (A) = \ n {B n ) and we 
want to determine P n so that B n is close to a diagonal operator at least for the entries 
with indices ranging between n — T n and n + r n where (t 11 )'^L 1 is a sequence of positive 
integers satisfying 

T n < Tn+i for n G N*, (4.2) 
t„ > oo, (4.3) 

n— >oo 

n - 2r„ ► oo. (4.4) 

n— *-oo 

We remark that in the proof of Theorem 1 1.1 1 we take r n = |_| n _| ' wnerc L S J := max{fc G 
Z : fc < s} means the integer part of s. 

Further on x G C 1 (M) is a fixed function satisfying 

0<X< 1, 

X(s) = 1 for s G [-1,1], 
X(*)=0forag [-2,2]. 

Then we write the decomposition 

Oi(j) = «n,/(j) + 5 n ,i(j) (4.5) 

with 

an,/0'):=«Ki)x(^), (4-6) 
a nJ (j):= ai (j)(l- X )(^-^ (4.7) 

and the corresponding decomposition 

A' = A n + A n , (4.8) 
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where 

A n = {S l a n ,i(A) + a n< i(A)S 1 *), (4.9) 

l<l<m 

A n = Yl {S l a n ,i(A) + a n< i(A)S 1 *). (4.10) 

l<l<m 

Using (|4.8I) we rewrite (|4.1[) in the form 

B n = c- iP " Ac iP " =D + A n + An + [D, iP n ] + W n , (4.11) 

where W n is considered as a lower order error. However due to 

n - 2r„ < j < n + 2r„ => a n> i(j) = (4.12) 

it is easy to see that for n large enough we have (D + A n )e n = De n = d(n)e n , i.e., d(n) 
is an eigenvalue of D + A n . Then in Section 14.21 Lemma 14.21 we show that d(n) is the 
n-th eigenvalue of D + A n provided n is large enough and the entries are sufficiently 
regular. Next we choose P n satisfying the commutator equation 

A n + [D,iP n ] = 0, (4.13) 

hence the expression (14. 1 1[) takes the form 

B n = e~ ip " Ae ip " = D + A n + W n 

and using X n (A) = \ n (B n ), d(n) — \ n {D + A n ) with the min-max principle we obtain 

\X n [A)-d{n)\ = \\ n (B n ) - X n (D + A n )\ < \\B n - (D + A n )\\ (4.14) 

for n > no- 

Lemma 4.1. Let P n be a finite rank self-adjoint operator satisfying A n = i[P n , D]. If 
A' = A n + A n and B n = e^ iP "Ae iP ", then 

\\B n -(D + A n ) || < ||[P„,i„]|| + i||[P„,A„]||. (4.15) 

All our results will follow from suitable estimates of the right-hand side of (|4.15l) . 
i.e., estimates of norms of commutators. A general estimate is stated in Section [5] and 
the norms of commutators from the right-hand side of (|4.15[) are estimated in Section [5] 

Proof. We introduce 

B n := e- iP "i„e iP " - i„ (4.16) 

and we observe that 

B n =J o Ys ( e " sP "^ eisP ") ds = l e~ isP "i[i„,P„]e isP "d s . 
Since for s G R the operators e lsP " are unitary, ||e IsPn || = 1 and we find 

\\B n \\ < ||[I nj P„]||. (4.17) 

Next for s £ R we introduce 

G n (s) := e _isP " (D + i[sP„, I?])e isPn — D (4.18) 
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and we observe that 

^G n {s) =e~ isP "(i[D + i[sP n ,D},P n } +i[P n ,D])e isP " 
= c- isP "s[[D 7 P n },P n ]e isP ". 

Therefore 

G n (l)= / e- isP - S [[D,P n ],P n }e isP "ds 
Jo 

and we can estimate 

||G„(1)|| < J 8\\[[D t P n ],P n ]\\ds = ±\\[[D,P n ],P n ]\\. (4.19) 

However 

B n = e- iP " (D + A n )e iP » + e - iP "i„e iP " = (G„(l) + D) + (B n + A n ), 

hence 

\\B n - (D + A n )\\ = \\B n + G n (l)\\ < \\B n \\ + \\G n (l)\\. (4.20) 

To complete the proof it remains to estimate the right-hand side of (I4.20j) using 1)4.170 - 

(I4T91) . " □ 

4.2. Equality d(n) — \ n (D + A n ). In this section we give sufficient conditions to 
ensure the equality d{n) = X n (D + A n ) used in estimate ()4.14|1 . 

We consider a sequence of positive integers (T n )^ =1 satisfying (|4.2|) - ()4.4|l and n > 0. 
We assume that the inequalities 

d(n) < d(n+l), (4.21) 

4mmax{|aj(n)|, \ai(n + m — t„)\} < d(n) — d(n + m — r n ) (4.22) 

hold for n > tiq and I — 1, ... ,m. 

Lemma 4.2. Assume that (|4.210 . 1)4.220 hold for n > no. If A n is defined by means 
of a n> i and \ as in Section \4- 1\ then there is hq G N such that 

d(n) = X n (D + An) for n>h . (4.23) 

Proof. We introduce 

OnU)~ J2 (\<lU)\ + ~ ( 4 - 24 ) 
l<l<m 

and observe that Lemma 13.21 allows us to estimate 

d- < \ n (D + A n )<d+ (4.25) 

with 

d- = inf Ki)-a(i)}, (4.26) 

d+=sup{d(j) + a(j)}. (4.27) 

i<« 

First step. We claim that 

j > n + r n a n (j) < d(J) - d{n) (4.28) 
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holds for n > hq + m. Indeed, replacing n by j — i in (|4.22p we obtain 

j - i>n Q =>■ 4m\ai(j - i)\ < d(j - i) - d(j -i + m- Tj-i), (4.29) 

hence for < i < m, j > n + r„ we have j — i + m — Tj_i > j — r„ > n and applying 
(|43T|) we find 

j > n + r n =+ d(j — i) — d(j — i + m — Tj-i) < d(j) — d(n) (4.30) 

for < i < m, n > n + m. Combining \d n ^(j)\ < \ai(j)\ with (|4.29[) and (14.301) we 
obtain ([Qg| . 

Second step. We claim that 

n + m < j < n - r„ => a n (j) < d(n) - d(j). (4-31) 

Indeed, replacing n by j — i + t„ in (I4.22[) with j > no + i we find 

4m|a ; (j - z) < d(j - i + Tj- i+Tn ) - d(j - i + m), (4.32) 

hence for < i < m, j < n — r n we have j — i + Tj-i +Tn < j + t„ < n and applying 
(|4~2T|) we find 

no + m < j < n — r„ => d(j — i + Tj-i +Tri ) — d(j — i + m) < d(n) — d(j) (4.33) 

for < i < m. As before, P~3"Tj) follows from (|4~3"2l and (|333")) . 

Next we observe that by definition n — t„ < j < n + t„ =+ a n (j) = 0, hence 

n +m<j <n a n (j) < d{n) - d(j) =+ d{j) + a n {j) < d{n), (4.34) 
n + m <n < j =*> d n (j) < d{j) - d{n) =+ d{n) < d(j) - a n (j). (4.35) 

Therefore one can choose no large enough to ensure d~ = din) = d+ for n > no and 
(|4~231) follows from (j4T231) . □ 

5. A GENERAL ESTIMATE 

5.1. Statement. We fix a sequence of positive integers (t„)^ 1 satisfying (|4.2p - (l4.4[) 

and for s > we denote 

a s (j) := max \ai(j + i)\, (5.1) 

l</<m 
|i|<s 

a n , s (j) ■= + max |(Ao,)0 + 01, (5-2) 

7V, K/<m 



|i|<s 



7,(j) : = , min (5.3) 

i|<s+l 

7.(j) :=max|(A 2 d)(j+z)|, (5.4) 

where (A 2 d)(n) = (Ad)(n + 1) - (Ad)(n) = d(n + 2)- 2d(n + 1) + d(n). 

Theorem 5.1 (general estimate). Let A be defined by (|1.1[) - (|1.6[) . Let (t„)^L 1; a s , 
7 S , a,i iS; % &e as above. If (I4.23|) ZioZds and 

/ -\ ~ /.N,a4mO) , - f .,a 4m (j) 2 

PnW ■= 2a n ,4mU) TV + m 75m(j) (5.5) 

74m 0) 74m 0) 
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then there is n\ such that for n > rt\ one has the estimate 

\X n (A) - d(n)\ < 15m 3 sup p n (n + i). (5.6) 

|'i|<2r„+4m 

Proof. This general estimate will be proved in Section [51 □ 
5.2. Application. We check that Theorem IBTTI implies 

Theorem 5.2. Let A be defined by Ql.l[) - (|1.6[) . Assume that there exist C > 0, 5 £ R, 
K > satisfying S + k > and «o £ N suc/i i/iai 

c -i n 5+K-i < (Ad)( n ) < Cn 5 ^" 1 /or n > n , (5.7) 

(A 2 d)(n) = 0(n 5+,i - 2 ) as n -»• oo, (5.8) 

<2/(n) = O^ 15 ) as n — > oo. Z = 1, . . . , m. (5-9) 

(Aaj)(n) = O^ 5-1 ) as n -> oo, 1 = l,...,m. (5.10) 

TTien one has the estimate 

X n (A) = d(n) + 0(n s ~ K ) as n -» oo. (5.11) 

Proof. Due to (|5.7[) and (|5.8[) there exist Co > cq > and ni e N satisfying 

c n S+K < d{n) < C Q n S+K for n > m. (5.12) 

In order to ensure (I4.23[) we will check that the assumptions of Lemma 14.21 hold if 
T n = n + m — [ n£ oJ where £o > is fixed sufficiently small. For this purpose we 
introduce c\ := Co — Co(so) S+K and estimate 

d(n) - d{lne \) > c n s+K - C (ne ) s+K = Cl n s+K for n > n x . (5.13) 

Let eo > be small enough to ensure c\ > 0. Then it is clear that (I4.22j) follows from 
(|5TT5j) and (HHJ). 

Thus all assumptions of Lemma l4.2l hold and it remains to apply Theorem l5.ll Using 
(|5.7I) - (|5.10I) we can find a constant C\ such that 

inf -f4m(n + i)>C^n S+K -\ (5.14) 

|i|<2-r„+4m 

sup 75m(n + i) < dn 5+K - 2 7 (5.15) 

|i|<2T„+4m 

sup o?4 m (n + i) < Girt" 5 , (5.16) 

|i|<2r„+4m 

sup <5 n ,4m (n + i) < Cirt" 5-1 (5-17) 

i|<2r„+4m 

with 7s, 7 S , a s , a rhSl given by (15.1l) - (|5.4p . Therefore, 

Wi) ~ C > V "^ = C ^"' <5 ' 18) 

hold when \ j — n\ < 2r„ + 4m and it is clear that 

sup p n {n + i) = 0{n s ~ K ) (5.20) 

|i| <2r„+4m 
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if pnij) is given by (1531) . We conclude that (|5.1ip follows from ([53)1 and (|5.20|) . □ 

5.3. Proof of Theorem 11.11 

Proof. Let <5 = 5 + n. Then the assumptions (HI) and (H2) imply (j5~7| - (j5~10|) . Con- 
sequently Theorem 11.11 follows from Theorem 15.21 □ 

5.4. Other applications of the general estimate. In this section we consider 
d(n) r~j uj(n) where the function lo : (0, oo) — > (0, oo) is one of a special type of functions 
described below. 

5.4.1. We fix k > 0, k' G R and assume 

\ai(n)\ + n\Aai(n)\ = 0(ri- K (lnn)- K ' u(n)), (5.21) 

d(n) =uj(n) (l + 0(n- 2 )) (5.22) 

where 

w(A) =c A' 5t, (lnA)' 5 o (5.23) 

holds with some Co > 0, 5q > 0, <5q G R. We observe that the derivatives satisfy 
o;( fc )(A) ~ CfcA _A: w(A) as A — > oo. Since Aw(n) = uj'(n + r n ) holds with some r n G [0, 1] 
and u/(n + r„) ~ o/(n) as n — > oo we easily deduce 

7s(n) ~ ^ n _1 a;(n) (n — > oo), (5-24) 
7 ^(n) = 0{n- 2 u{n)) (n -> oo). (5.25) 

Using t„ = L n /4J we find that Theorem 15. II gives the estimate 

X n (A) = d(n) (l + 0(n~ 2K (lnn)- 2K ')) . (5.26) 
It is easy to see that (|5.26[) still holds when k = and n' > 0. 

5.4.2. We assume that (|5.21[) holds with some k > 0, k' G E, 

<2(n) = w(n) (l + 0(n- 2 (lnA)" 1 )) (5.27) 

where 

w(A) = c (lnA) 5 o (5.28) 
holds with some Co > 0, (5q > 0. Then computing the derivatives of lo we find 

7s(n) ~ <5on _1 (lnn)~ 1 a;(n) (n — >■ oo), (5.29) 
7,(n) = 0(?i" 2 (lnn)- 1 o;(n)) (n4 oo); (5.30) 

and using r„ = L"/4J in Theorem 1 5. II we obtain the estimate 

X n (A) = d(n) (l + 0(n- 2K (lnn) 1 - K ')) ■ (5-31) 

It is easy to see that (|5.31[) still holds when n = and k' > 1. 
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5.4.3. We assume that ft > 0, < 6> < 1 and 

|aj(n)| + n 1_e |Aai(n)| = 0(n- K uj(n)), (5.32) 

d(n) = u(n) (1 + 0{n 9 - 2 )) (5.33) 

where 

W (A) = c Q X 5 ° e cA9 (5.34) 

holds for some c > 0, c > 0, S <E R. The derivatives satisfy o/ fc )(A) ~ c k \^ k ^^ 9 ^uj(\) 
as A — > oo and we deduce 

ls (n) ~ c6n- {1 - e) uj{n) (n ->• oo), (5.35) 

7 ^( n ) = 0(n" 2(1 " e) w(n)) (n oo). (5.36) 

Let t„ = L"- 1_e £oJ with £o > small enough. Then d(n) — d(n + m— r n ) > d(n)/2 for 
n > ri\ and Theoren l5 . 1 1 ensures 

X n (A) = d(n) (1 + 0(n- 2K )) . (5.37) 

6. Proof of Theorem 15.11 

6.1. Estimates of commutators. For n 6 N*, 2 e Z we consider a„ ; : Z — > K 
satisfying 

Pn-l(j) = Pn,l(j) and a n-l(j) = On,l{j)> ( 6A ) 

Pn,l(J) = a n,lU) = when 3 ^ °> ( 6 - 2 ) 

Pn,l(j) = 07i,i(j) = when \l\ > m, (6.3) 
where m € N* is fixed. We assume moreover 

Pn,lU) = wnen li _ n \ ^ 2t «> ( 6 - 4 ) 
where (T n )^ =1 is as before and consider finite rank self-adjoint operators 

P n -= ]T (S l p n ,i(A)+p n ,i(A)S 1 *), (6.5) 

l<l<m 

R n := i[P n , A n ] = i(F„A„ - A n P n ), (6.6) 

where A n is given by (I4.9p . For s > 0, j G Z define 

a«, s (j) := max \a n .i(j + i)\, (6.7) 
|?|<m 
|i|<« 

Pn.s(j) ■= max \p n .i(j (6.8) 

|(| <rn 
\i\<s 

«n,s CO : = max I ( Aa «,' ) C? + 1 > ( 6 - 9 ) 
|i|<m 

|i|<« 

^(^^maxKAp^Oy + i)!. (6.10) 

|t|<m 
\i\<s 
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Lemma 6.1. Let A n , P n , R n , a„ iS , l3 n , s , ct' n s , j3' n s be as above and 

Pn,s(j) ■= an,s(j)P'n,sU) + a n, a U)Pn,s(J)' ( 6 - U ) 

Then one has 

\\Rn\\<lOm 3 sup PnAm(n + i). (6.12) 

|i|<2T„+4m 

Proof. We can express R n as 

Rn = r n , (A)+ J2 (S k r n , k (^)+rnA^S k *) , (6.13) 

KK2m 

where for k > we have 

r n ,k(j) = (ej+k,R n e 3 ) = i(e J+fc , P„A„ e^) - i{e j+k ,A n P n Cj). (6.14) 

For s£l we write s+ = max{s, 0} and s_ = (— s) + . Then 

(e i; P n Cj) = p n ,i-j(i j)+), (6.15) 

(e l ,A n e j ) = a n j-j(i - (i - j)+), (6.16) 

and using (|6.15[) . (|6.16[) in 

i 

£j+k,A n ej + k-i)(e 

I 

we find 

r„,fe(i)=i ^2 r ntk ,i(j) (6.17) 

l<|i|<m 

with 

r n ,k,i(j) =Pn,k-l(j + k - (k - l) + )a rii i(j - Z_) 

- a n ,l{j + k- l+)Pn,k-i(j -(k- I)-). 
Moreover p n .k-i => \k — l\ < m and we claim that 

r n ,k,iU) ^ U - »l < 2r » + 2m. (6.18) 
Indeed, it suffices to use < fc < 2m and \k — l\ < m in 

Pn,k~i(j - (k - I)-) ^ =>• n - 2t„ < j - (k - I)- < n + 2t„, 
Pn,k-i{j + k— (k — l)+) ^ =>• n - 2r„ < j + k - (k - l) + < n + 2r„. 
Then reasoning as in Section [2] we can estimate 

\\Rn\\ < ||r„,o(A)||+ ]T (2|k„, fc (A)|| 2 +2||r„, fe (A-fc)|| 2 ) 1/2 (6.19) 

Ki<2m 

and the right-hand side of (|6 . can be estimated by 

(4m + 1) sup max \r n ,k(j — (6.20) 

j~>l 0<i</c<2m 

However for i, i' £ Z such that i' < i we have the expression 

a n ,i(j + i)-a„,i(j + i')= J2 (Aa„.;)(.?+A (6.21) 

i'<j'<i-l 



ASYMPTOTICS OF LARGE EIGENVALUES FOR. A CLASS OF BAND MATRICES 



15 



and using \k — Z+ + Z_| = \k — l\ < m we obtain the estimate 

\a n ,i(j +k-l + ) - a n ,i(j < ma' ns (j) (6.22) 
with s = max{|Z|, \k — l+\} < max{m, k} < 2m. Thus (|6.22p holds with s = 2m and 

\p n ,k-i{j + k-{k- 0+) - Pn,k-l{j ~(k- l)-)\ < mp' n;2m (j) (6.23) 
follows similarly. Since r n ^{j) = r' n k (j) + r" k (j) holds with 

r' n ,kU) = (Pn,fc-i(j + k-(k- 0+) -Pn,k-l{j - (k - Z)_))a n ,,(j - *-), ( 6 - 24 ) 

< fc (i) =p n ,*-iCj - ( fc - 0-)(on,i(j - *-) - + *- *+))> (6-25) 

we obtain |r nj fc(j)| < m p n ^ m { j) and ||i? n || can be estimated by 

2m 2 (4m+l) sup max \p nt 2 m (j — (6.26) 

|j-n|<2r n +4m 0<i<2m 

To complete the proof of (|6.12p it remains to use p n ,2m{j — i) < Pn,2m+|i|(i)- 1— ' 
6.2. Proof of Theorem [HU 

Proof. Let no be as in (|4.21|) and for j > 1 — / > no denote 

rfRi):=d(i + /)-d(i). (6.27) 
Consider A„, P n given by (TO)) . (TCTu]) with a„ : ; as in (fH))) and 

. an,i{j) 

Then R n := [iP„, D] coincides with A n due to 

i(e i+ i, R n ej) = i(e j+ i, P n Dej) - i(De j+ i, P n ej) 
= i(d(j)-d(j + l))(e j+l , P n e 3 ) 
= a n ,i(j) = (e j+ i,A„ej) for / > 0. 

Thus due to (|4~T4l and (l4"T5)l 

\X n (A)-d(n)\ < ||L4„,P„]|| + i|P„,P„]|| for n>n . (6.29) 



Pn.i(j) = ( e j+h Pne 3 ) =i^777- for Z = 1, . . . , m. (6.28) 



We consider a„ iS , /3„ )S , a4, s , (3' nyS given by (|6. T[) - (|6 . 10|) and in order to apply Lemma ISTTl 
we will check that 

PnAmti) < Pn(j) for n > n + 5m (6.30) 
holds with p n given by (|5.5|) and by (|6.11l) . Indeed, it is clear that 

a n ,s(j) < a s (j), (6.31) 

/3„, s (i) < AO') : = ^44 for i > ™o + m + s (6.32) 

hold with 7 S , a s given by (|5.3[) . (|5.1[) . Then we observe that the function x considered 
in Section 14.11 can be chosen such that 

llx'lloo : = suplx'WI <2, 

where denotes the derivative of x an d consequently 

|X((? + 1 - n)/r n ) - X ((j - n)/r n )| < 2/r„, 
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hence taking a„ iS as indicated in (|5.2p we obtain 

<«(?*) < fi»,-C?) ( 6 - 33 ) 
due to (|4.6I) . Finally we need to estimate f}' n s (j) and for this purpose we first observe 
that 

j > n Q & I > d{(j) > d(j + 1) - d(j) > 0, 
3 > n Q - I k I < => -d[(j) > d(j) - d(j - 1) > 
hold due to (|4.21l) . Further on we assume j > no + m. Thus for |/| < m we have 

mm{K(j)|, l4(i + l)|}>7o(j) (6.34) 
and writing (Ad;)(j) = E 4e i(i) A 2 cZ(j + i) where Z(Z) = [0, / - 1] n Z when / > and 
X(Z) = [I + 1, 0] n Z when Z < 0, we obtain 

< mm- 

Hence, 

1 1 (AdjKj) 



dKi + i) W 

However by definition (|6.28p we have 



d{(7 + l)d{Ci) 



< 



a I -m ^ l Aa »j(j)l ,, 



K(i + i)l 



A 



7o(j) 2 ' 
1 



(6.35) 
(6.36) 

(6.37) 



hence using \l\ < m, (|6.34l) and f|6 .36(1 to estimate the right-hand side of (16.37[) we 
obtain 

< BnAj) ^rr + «.(?) (6.38) 



js(j) is{jY 

for n>n + m + s. Now it is clear that (1630)) follows from ([QTj) - (jQ3"]) and ([6351 . 
Then using (|6.30[) in Lemma T3. 2 1 we obtain 

|[A„,P„]| < IOto 3 sup p n (n + i) for n> n\. (6.39) 

|i|<2r„+4m 

To complete the proof it remains to observe that one can use a s and a n . s as above if a n j 
replaces a n> i, hence the norm |[A„, P n ]\\ can be estimated in a similar manner. Thus 
we can conclude that the right-hand side of (|6.29p can be estimated by the right-hand 
side of (EH). □ 



Acknowledgments. The second named author's research was partially supported by 
TODEQ MTKD-CT-2005-030042. 

References 

[1] A. Boutet de Monvel, S. Naboko, and L. O. Silva, The asymptotic behavior of eigenvalues of a 
modified Jaynes-Cummings model, Asymptot. Anal. 47 (2006), no. 3-4, 291—315. 

[2] A. Boutet de Monvel and L. Zielinski, Explicit error estimates for eigenvalues of some unbounded 
Jacobi matrices, Spectral Theory, Mathematical System Theory, Evolution Equations, Differential 
and Difference Equations: 1WOTA10, Oper. Theory Adv. Appl., vol. 221, Birkhauser Verlag, 
Basel, 2012, pp. 187-215. 

[3] H. Behncke and F. O. Nyamwala, Spectral theory of difference operators with almost constant 
coefficients, J. Difference Equ. Appl. 17 (2011), no. 5, 677-695. 

[4] , Spectral theory of difference operators with almost constant coefficients II, J. Difference 

Equ. Appl. 17 (2011), no. 5, 821-829. 



ASYMPTOTICS OF LARGE EIGENVALUES FOR A CLASS OF BAND MATRICES 



17 



[5] J. Janas and S. Naboko, Infinite Jacobi matrices with unbounded entries: asymptotics of eigen- 
values and the transformation operator approach, SIAM J. Math. Anal. 36 (2004), no. 2, 643-658. 
[6] M. Malejki, Asymptotics of large eigenvalues for some discrete unbounded Jacobi matrices, Linear 

Algebra Appl. 431 (2009), no. 10, 1952-1970. 
[7] , Asymptotic behaviour and approximation of eigenvalues for unbounded block Jacobi ma- 
trices, Opuscula Math. 30 (2010), no. 3, 311-330. 
[8] M. Reed and B. Simon, Methods of modern mathematical physics. II. Fourier analysis, self- 

adjointness, Academic Press [Harcourt Brace Jovanovich Publishers], New York, 1975. 
[9] M. A. Shubin, Pseudodifference operators and their Green function, Izv. Akad. Nauk SSSR Scr. 
Mat. 49 (1985), no. 3, 652-671. 
[10] G. Teschl, Jacobi operators and completely integrable nonlinear lattices, Mathematical Surveys 

and Monographs, vol. 72, American Mathematical Society, Providence, RI, 2000. 
[11] E. A. Tur, Jaynes-Cummings model: solution without rotating wave approximation, Optics and 

Spectroscopy 89 (2000), no. 4, 574-588. 
[12] H. Volkmer, Error estimates for Rayleigh-Ritz approximations of eigenvalues and eigenf unctions 

of the Mathieu and spheroidal wave equation, Constr. Approx. 20 (2004), no. 1, 39-54. 
[13] L. Zielinski, Eigenvalue asymptotics for a class of Jacobi matrices, Hot topics in operator theory, 
Conference Proceedings, Timi§oara, June 29-July 4, 2006, Theta Ser. Adv. Math., vol. 9, Theta, 
Bucharest, 2008, pp. 217-229. 

1 Institut de Mathematiques de Jussieu, Universite Paris Diderot Paris 7, 175 rue du 
Chevaleret, 75013 Paris, France, E-mail: aboutet@math.jussieu.fr 

2 Instytut Matematyczny PAN, ul. Sw. Tomasza 30, 31-027 Krakow, Poland, E-mail: 
najanas@cyf-kr.edu.pl 

3 LMPA, Universite du Littoral, 50 Rue F. Buisson, B.P. 699, 62228 Calais, France, E-mail: 
Lech.Zielinski@lmpa.univ-littoral.fr 



